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A recent paper by Martin and Schwarz |y argues that the "standard infla- 
tionary result" has been finally proven. The result itself is formulated as: "the 
closer the inflationary epoch is to the de Sitter space-time, the less important 
are large-scale gravitational waves in the CMBR today". Beginning from 
the basic equations of Grishchuk H , || the authors say [1] that Grishchuk's 
conclusion about approximate equality of metric amplitudes for gravitational 
waves and density perturbations "is wrong because the time evolution of the 
scalar metric perturbation through the (smooth) reheating transition was not 
calculated correctly". They reiterate a claim about "big amplification" of 
i ■ scalar perturbations (in contrast to gravitational waves) during reheating. 

(3JT), The authors say [1] that after appropriate correction they have recovered the 

"standard result" within Grishchuk's approach. It is shown in this Comment 
that the "big amplification" is a misinterpretation. There is no difference in 
jJJ ■ the evolution of long- wavelength metric perturbations for gravitational waves 

and density perturbations: they both stay approximately constant. The in- 
fluence of cosmological transitions on the evolution is none at all, as long as 
the wavelength of the perturbation is much larger than the Hubble radius. It 
is shown that from the approach of [2, 3] follow the conclusions of [2, 3] with- 
out change. Finally, it is argued that the "standard inflationary result" does 
not follow from the correct evolution and quantum normalization of density 
perturbations. 
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Cosmological gravitational waves and density perturbations could have been generated 
as a result of the superadiabatic (parametric) amplification of their zero-point quantum 
oscillations by strong variable gravitational field of the early Universe. The theoretical 
predictions regarding the amplitudes and spectral slopes of these primordial perturbations 
are of a considerable practical importance, especially in view of the forthcoming sensitive 
observations. 

One possibility is that the early Universe was driven by a scalar field (p, which has 
latter decayed into a normal matter. For this model, some theoretical studies of density 
perturbations have led to the what is known as the "standard inflationary result". The 
central inflationary formula for the generated matter density variation is usually written as 

(8p/p) f ~ (H 2 /<p) i 

where the left-hand-side is what we want to know at the end of the long-wavelength regime 
(second Hubble radius crossing), and the right-hand-side is what we need to evaluate at the 
beginning of the long- wavelength regime (first Hubble radius crossing). The same formula 
is often written as 

(5p/p) f ~ (V^/V^ 

where V((p) is the scalar field potential and V — 

To say the least, this formula looks very counter-intuitive. Imagine that the potential 
V(ip) has an inflection point, where V = 0, and the slowly-rolling scalar field passes at 
a certain moment of time through that point. Then, the final amplitude (5p/p)f of the 
mode, which entered the long-wavelength regime at that moment of time, is predicted to be 
infinite. Since the amplitude of the metric perturbation hf at the second crossing is of the 
same order of magnitude as the matter density variation (5p/p)f (both are dimensionless 
numbers), the metric (curvature) perturbation is also predicted to be infinite. Moreover, 
for the generation of a spectral interval with the almost Harrison- Zeldovich slope n = 1 
one needs an interval of the almost de Sitter expansion, that is, the potential V((p) should 
be very flat, V ~ 0, at some interval of <p. Then, all the amplitudes of the generated 
density perturbations in this interval of spectrum are predicted, according to the "standard" 
formula, to be arbitrarily large. In addition, all these predictions are in a sharp conflict with 
the generation of gravitational waves, whose final amplitudes are finite and small, even if 
the discussed interval of the early expansion was of the de Sitter type. 

There is not any obvious physical reason for the appearance of this divergent result in a 
space-time (gravitational pump field) with finite and relatively small background curvature 
(the Hubble parameter Hi is probably 5 orders of magnitude smaller than the Planckian 
value) which sets the numerical level of the nominator in the "standard" formula. This 
result is also in conflict with the finite "temperature" || of the de Sitter space-time, which is 
determined by its constant Hubble parameter H. Whatever "particles" are being produced, 
they are supposed to have finite and small energy in every frequency interval. However, 
according to the "standard" formula, even a short interval of the de Sitter evolution is capable 
of generating an arbitrarily large amount of "particles" - density (curvature) perturbations. 

The results of calculations of ref. [2, 3] have shown that the "standard" formula does 
not follow from the proper evolution and quantum normalization of density perturbations. 
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Although the calculations are quite involved technically and include elements of quantum 
mechanics, the bottom line is simple and can be summarised as follows. Cosmological den- 
sity perturbations are not perturbations in the matter component only, be it a scalar field 
or a perfect fluid; they necessarily include the metric (gravitational field) perturbations, and 
do not exist without them. We are mostly interested in the evolution of the metric pertur- 
bations, because they are the carriers of density perturbations - the scalar field, together 
with its perturbations, is supposed to be converted into other forms of matter anyway. The 
evolution of the metric perturbations is trivial: the dominant solution is a constant, as long 
as the wavelength is much larger than the Hubble radius, and indepenedently of what kind 
of transformations of the matter content are taking place during all the time of expansion 
from the beginning of the long-wavelength regime and up to its end. So, the final amplitude 
of the metric perturbation is practically the same as the initial one, hf ~ hi. To predict the 
final amplitude we need to know the initial amplitude. Classically, the initial amplitude is 
in our hands and is determined by our choice of initial conditions. Quantum-mechanically, 
the minimal initial amplitude hi is determined by quantum fluctuations. To find the nu- 
merical value of hi we certainly cannot rely (in contrast to what is often being done) on 
the quantization of a free scalar field, which is not accompanied by metric perturbations 
at all. We need to quantize the "fundamental" field, which fully represents cosmological 
density perturbations and includes metric perturbations. The dynamical equation for this 
"fundamental" field is very similar (as well as the Lagrangian, Hamiltonian and other struc- 
tures) to the dynamical equation for cosmological gravitational waves, and they are exactly 
the same in case of any of the power-law expansions. The initial amplitude hi, determined 
by the quantum fluctuations for a given mode with wavelength A, is hi « lpi/Xi, where 
Xi = c/Hi and Hi is the Hubble parameter at the first crossing. This finite and small 
number effectively translates into hf and (5p/p)f, that is, (5p/p)f ~ hf ss hi. Since the 
initial amplitude for gravitational waves (derived from its own quantum considerations) is 
also approximately equal to hi, and since it is being transmitted without change to the end 
of the long-wavelength regime, the final amplitudes of gravitational waves and density per- 
turbations are also approximately equal. It follows from the calculations of ref. [2, 3] that 
in the correct version of the "standard" formula the right-hand-side of this formula must 

be multiplied by the dimensionless factor (\J —H / H 2 )i. This factor cancels out the zero in 
the denominator and makes the amplitudes of density perturbations finite and of the same 
order of magnitude as the amplitudes of gravitational waves. 

Martin and Schwarz [1] (as well as some of their predecessors) disagree with the outlined 
above and used in [2, 3] law of evolution for the metric perturbations associated with den- 
sity perturbations. They claim that the long-wavelength scalar metric perturbations have 
experienced, in contrast to the metric perturbations representing gravitational waves, "big 
amplification" during "reheating", that is, during a short interval of time when the scalar 
field decayed into the radiation-dominated matter. In other words, they try to explain the 
huge number for hf suggested by the "standard" formula by the evolution of a small initial 
number which was determined by quantum fluctuations. This seems to be the only way to 
make sense of the "standard" formula. Otherwise, if you agree with the law hf ~ hi, the 
huge hi must be postulated from the very beginning, and this is certainly not the amplitude 
dictated by the zero-point quantum oscillations of a quantized field. The evolution of the 
long- wavelength metric perturbations is the main point of the disagreement, and we will 
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address it below. 

Let us formulate the problem more precisely. We consider classical solutions to the 
perturbed Einstein equations. The equation of state of the background matter may vary 
as a function of time. A known example is the decoupling transition from the radiation- 
dominated era to the matter-dominated era. The decoupling transition is not an issue of 
much disagreement, so we will concentrate on the reheating transition from the era governed 
by the scalar field to the era governed by the radiation. Theoretically, the transition could 
be either smooth, so that the scale factor, its first time derivative, and its second time 
derivative are continuous; or sharp, so that the second time derivative of the scale factor 
(and the first time derivative of the Hubble parameter) experience a finite jump. Other 
quantities of the background solution are either nonimportant or calculable. We consider 
density perturbations (scalar eigen-functions in Lifshitz's classification) and gravitational 
waves (tensor eigen-functions in Lifshitz's classification). 

The question is: What can be stated about the evolution of metric perturbations in 
the long-wavelength regime ? The provisional answer to this question, even if not for all 
conceivable cases, is effectively known since the pioneering paper of Lifshitz || : the dominant 
solution is a constant. The perturbation in the matter energy density is slowly growing, all 
the time from the first Hubble radius crossing to the second Hubble radius crossing. But 
the perturbation in the metric tensor remains constant, both, for density perturbations and 
for gravitational waves. So, if we have started deeply in the very early Universe with the 
perturbed metric amplitude at the level, say, 10 -4 , we will get the same number 10 -4 at 
the end of the long-wavelength regime. The end of this regime is today for wavelengths 
of the order of today's Hubble radius. And the initial metric amplitude, say, 10 -20 also 
translates into 10~ 20 today, both for gravitational waves and density perturbations, and 
nothing dramatic happens on the way from one crossing to another. The perturbation in 
the matter energy density (in case of density perturbations) will grow, but only up to the 
level 10~ 4 or, correspondingly, 10 -20 . 

This answer requires a few clarifications. Which metric perturbations are we work- 
ing with ? After a few attempts one usually realizes that the class of the battle-tested 
synchronous (freely falling) coordinate systems is a very convenient choice for practical cal- 
culations. There exists some remaining freedom in the choice of a synchronous coordinate 
system. This freedom is an advantage rather than disadvantage, because it allows us, if we 
want, to adjust the coordinate system in such a manner that by the time of reaching the 
matter-dominated stage it will become a comoving coordinate system in addition to already 
being synchronous - the most favorite choice. However, we will mostly work with a single 
function, denoted below as p,/a, which completely describes the metric perturbations, and 
which is independent of the remaining freedom. It is for this function we gave the estimate 
of hi and hf ~ hi. In general, there are two essential metric components (two polarisation 
states) for each type of perturbations. These components are independent functions in case 
of gravitational waves. But for the density perturbations, that we actually consider, one po- 
larisation component (longitudinal-longitudinal) is not independent and can be expressed in 
terms of the another (scalar) component. From a given complete solution one can calculate 
any of the so-called "gauge-invariant" variables, and we will do this calculation below for 
some of them. 

As was already mentioned, the dominant ( "growing" ) long- wavelength solution for each 
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independent metric component, both in case of gravitational waves and density perturba- 
tions, is a constant C. For simple numerical evaluations, it is convenient to operate with 
a single number - the characteristic amplitude - which is the square root of sum of squares 
of the two polarization components. In practice, the characteristic amplitude is again C. 
For the perturbations with wavelengths of the order of today's Hubble radius, this number 
C gives also a rough numerical estimate |J, [0 of the produced quadrupole anisotropy in 
the CMB radiation: 5T/T = C. And, finally, if this constant C is approximately the same 
number for gravitational waves and for density perturbations, the anisotropics 5T/T induced 
by gravitational waves and by density perturbations will also be approximately equal. Why 
these numbers C should be approximately equal for primordial density perturbations and 
for primordial gravitational waves is determined by the initial conditions. This is the only 
one place where quantum considerations relate to our discussion. We have mentioned this 
point above and will return to it later, but it is not in the focus of our attention to the 
claimed dramatic "influence of transitions" in the case of density perturbations. 

In a sense, we know in advance the correct answer to the question of evolution of long- 
wavelength metric perturbations: the dominant solution is a constant. This expected answer 
was essentially rederived in great detail in || and has led to the conclusion that the "stan- 
dard" result does not follow from the correct evolution and quantum normalization of density 
perturbations. The work [2, 3] had been several times critisized and "refuted" by members of 
inflationary community. Martin and Schwarz |l|] make an analysis of two previous dedicated 
papers |§, @, and we will need their summary. The authors of [0 note that "according to 
Deruelle and Mukhanov, Grishchuk made two mistakes: he took wrong joining conditions 
and he used the wrong equation of state. ..at the reheating transition". However, Martin and 
Schwarz conclude that "both his joining conditions at the reheating transition and the 
equation of state after reheating have been used correctly". Martin and Schwarz [[I]] recall 
that according to Caldwell, Grishchuk's conclusion occurs because the long-wavelength limit 
n 2 —>■ "has not been taken consistently". However, they conclude, "the limit had been 
taken properly" by Grishchuk. So, we can now safely put aside other papers and concentrate 
on the "final proof Q of the statement that, nevertheless, Grishchuk's conclusion is wrong 
because of the overlooked "big amplification" of the metric perturbations during reheating. 

We need some equations and relationships. Most of them are reproduced in |3J but it is 
better to have them, so to say, from the first hands. 

We write the perturbed metric in the form 

ds 2 = -a 2 (?7)[d?7 2 - (5ij + h ij )dx i dx j ]. 

We will often use two new functions of the scale factor a(rj): a = a' /a and 7 = 1 — a' jo? = 
1 + (a/ a')'. In terms of t time, cdt = a(r])dr], the function 7 is 7 = —H/H 2 where H is 
the Hubble parameter H = a/a = ca/a. As a consequence of the background Einstein 
equations, we have 1 +po/eo = (2/3)7. F° r models under consideration, the instantaneous 
equation of state may vary, smoothly or sharply, from p Q = — e (7 = 0) to p = e o/3 (7 = 2) 
and later to p = (7 = 3/2). 

In case of density perturbations, we write a spatial Fourier component of the metric 
perturbations in the form 

hj = h(rj)SijQ + hi{rj)—Q,ij 
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where Q = e m x or e _mx . The function h(rj) represents the scalar polarisation state, whereas 
the hi(rj) represents the longitudinal-longitudinal polarisation state. 

It was shown [2, 3] that, for models governed by scalar fields and perfect fluids, the 
full set of perturbed Einstein equations can be reduced to a single second-order differential 
equation (master equation). To write this equation specifically for models governed by a 
scalar field with arbitrary scalar field potential, we introduce the function /x/a: 



^ =/! + -/!' 

a aj 

and the function /x: /x = y/jp,. Then, the mentioned master equation takes the form 



(1) 



fJ," + fJL 



n 



aV7 



0. 



(2) 



The function jx/a can be called the "residual-gauge-invariant" part of h(rj). This function 
plays the central role in our discussion. 

All the functions describing the perturbations - h(r]), hi(r)), and (fi(r)), where (pi(rj) is 
the scalar field perturbation, <p = <fo(v) + Vii^Q ~ can now be found from solutions to eq. 
(§). In particular, using (|l|) we find 



a 
a 



VO CI 



(3) 



where is an arbitrary integration constant reflecting the remaining coordinate freedom. 
Equation (0) has been presented, on purpose, in the form similar to the one of the 



previously explored equation for gravitational waves |TC| : 



//' + /x 



n 



a" 
a 



where, in case of gravitational waves, 



h - Vld 



(4) 



(5) 



and Qij = p^e 111 '*-, Pij5 1 ^ = 0, Pijn^ = for each of two independent polarisation tensors Pij. 
(Equations (0) and (f|) are reproduced as (3.24) and (2.16) in [1].) 

Equations and solutions for the function /x representing density perturbations can be 
found from equations and solutions for the function /x representing gravitational waves by 
the simple replacement a — > [2]. 

We are now in the position to find long- wavelength solutions to (^) and @. It is known 
lOH that as long as the wave stays "under the barrier", that is n 2 -C \a"/a\ 
solution to (W) is 



the approximate 



f di] 

/x = da + C 2 a j — 



(6) 



where C\ and C2 are arbitrary constants. Combining (|6|) with (^) and neglecting the decaying 
term with C2, we find the dominant solution for the gravity- wave characteristic amplitude 
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h 9 - w - = C v (7) 

In other words, independently of the sort of matter which drives the scale factor a(rj) of our 
model, the dominant solution is a constant during all the time that the wavelength is longer 
than the Hubble radius. (A certain distinction between the notions of the "barrier" and the 
Hubble radius "crossing" is not important for our discussion.) 

Similarly to the gravity- wave case, the long-wavelength solution to eq. @ is 0: 



fi = C ia + C 2 a [-£-, (8) 
and the dominant part is /i oc 0^7, that is, 



£ = C, (9) 

and hi = hf = C\. Again, the scale factor a{r\) and the function 7(77) change with time, 
and, in particular, the function 7(77) may be changing from very small values up to 7 = 2, 
but solution @ is still valid. Using @ and the fact that 07 = (a/a)', we find from (|3|): 

h(r)) = d + 

a 

As was explained above, the term with Cj is useful, but in any case this term is decaying and 
we neglect it for the purposes of our discussion. It can also be shown [[| that the component 
hi(rj) slowly varies but is small all the time, and can reach the numerical level G\ only at the 
end of the long-wavelength regime. So, in case of density perturbations, similarly to eq. (^) 
for gravitational waves, the dominant solution for the characteristic metric amplitude is a 
constant 

h d - p - = Q. (10) 

This is true independently of whether the scale factor 0(77) is driven by a scalar field or by a 
perfect fluid. Of course, the constancy of the metric amplitude at the radiation-dominated 
and matter-dominated stages is known for long time JO]] . 

The problem of the evolution admits also a complete solution in case of a sharp transition. 
Imagine that the function 7(77) is discontinuous at the transition point and takes the values 
7; and 7 r to the left and to the right of the transition point. The ji can be arbitrarily close 
to zero, and 7 r can be 2. It was shown [FJ] that eq. (0) requires the continuity of the function 
v, where 

v = 7 fiy (ii) 

a 

and, hence, the continuity of the function p,/a. The continuity of jx/a can be derived either 
directly from (0) or from inspection of (PP. Indeed, a jump of Ji/a would develop a 5 - 
function which could not be compensated by the step function 7 in the right-hand-side of 
(|TT|), and hence the continuity of v would be violated. Instead, eq. flTTp requires a finite jump 
in the first derivative of Ji/a. The value of this jump is determined by the jump of 7 and 
by the continuity of v ||12|| . Obviously, the continuity of Ji/a guarantees that the constant 
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G\ is being carried over any transition, including the reheating transition from po = —6q 
to po = eo/3, without any changes. The "influence of reheating", whether the transition is 
smooth or sharp, on the scalar metric perturbations is none at all. 

At this point it is necessary to say that the formula fi\ v=ril -o = fi\ v=Vl+ o printed at p. 7161 
of [2] is indeed an error (which could have caused trouble for Martin and Schwarz). [There is 



also another obvious error on the same line: J (2 + + 0) should read (2 + 0)/{l + /?).] 
In the joining condition, it was meant to be ft, not /i, but the error arose because of the 
change of notations in original versions of the manuscript. It is obvious that this error had 
not been used and had not influenced the calculations. It is clearly stated at the very end of 
p. 7164 which condition had actually been used: "From the continuous joining of h(rj). ..one 
derives...", which implies the continuous joining of the function fl/a, not \i/a. The full set of 
conditions is reproduced as eq. (5.11) in ref. [1], which, the authors agree [1], "are exactly 
the conditions that have been used" in ref. [2]. So, there was no mistake in calculations, 
caused by the mentioned error at p. 7161 of [2], contrary to what the authors of [1] suspect 
in their Sec. VI. Their eq. (6.5) is correct and the same as was actually used in [2], but it 
does not lead to any overlooked "big amplification" as we shall see below. 

We have shown that the evolution of scalar metric perturbations is as simple as the evo- 
lution of tensor metric perturbations. This evolution is at the basis of practical calculations 
in [2, 3]. There is no room for a "big amplification", which would pick up specially the 
metric amplitude of density perturbations and increase it, say, from 10~ 20 to 10 -4 , while 
leaving the metric amplitude of gravitational waves at the initial level 10~ 20 . We gave the 
complete solution to the problem of evolution and did not need any concept of the gauge- 
invariant formalism or ( conservation law. But in order to demonstrate that the claimed 
"big amplification" is a misinterpretation, we need to deal with those concepts. 

Martin and Schwarz |l| work with the gauge-invariant potentials $ and (. The important 
role is allocated to the "constancy of (" . The quantity ( is defined by their eq. (4.8): 

2 H" 1 ^' + $ 

C = - (12) 
S 3 1+w ' V ; 

where Ti = a and l + w = l +po/eo = (2/3)7. 

Since we have solved our perturbation problem, any quantity, including any gauge- 
invariant potential, can be calculated. It was shown [2b, 3] that the definition of $ (see 
(2.3) in [1]), plus available Einstein equations, produces the exact equality 

• = 5?°*!>' (13) 



(reproduced as (3.27) in [1]). It was also shown [2b, 3] that the definition (|T2"D, plus ( |i~3|) 
and, most importantly, plus the master equation (Q), produces the exact equality 

C = -i(£). (M) 

In other words, after properly taking into account all the equations, the quantity £ turns out 
to be simply the function fl/a (up to the factor —(1/2)) which we are working with. Equality 
( [14]) demonstrates that the ( conservation law is empty, in the sense that it does not give us 
anything new above what we have already derived in a much simpler way [2b, 3]. We know 
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already (see @) that the growing part of the function fl/a, in its lowest nonvanishing long- 
wavelength approximation, is a constant C\ - this fact immediately follows from equation 
(0). If one likes, one is free to work with the function (, facing the same problem that we 
discuss here. One should give physical interpretation to the function ( and realize that the 
"standard" formula predicts a huge value for this quantity at the second crossing, because 
Q ~ (5p/p)f. Since (f ~ Q, this huge value should have been postulated from the very 
beginning, from the first crossing. Obviously, the ( conservation law, which the authors [1] 
defend so strongly, does not allow any "big amplification" of ( on the way from one crossing 
to another. 

This appears to be the end of the story about constancy of (. However, we need to 
return to its beginning in view of the discussion in [1]. The fundamental equation often used 
in the literature is (see eq. (2.11) in [1]): 

$" + 2 (H - ^)$' + [n 2 + 2(W - H^)}$ = 0. (15) 

The original derivation and use of the "constancy of ( for superhorizon modes" was based 
on eq. (|15j) in which the term n 2 $ was neglected. Then, in terms of ( defined by ([12]), this 
truncated equation is equivalent to (' = from which the "conservation law" 

( = const = X 

was first derived (see for example [0]). It was shown [2b, 3] that the fundamental equa- 
tion (0) requires this constant X to be a strict zero, this constant X is not determined by 
initial coditions for the function fi(r]). The constant X arises because eq. fll~5|) is a third- 
order differential equation in terms of /j,(t)). In the decomposition over small parameter n 2 
this constant X would stand in front of n~ 2 : ((rj) = X/2n 2 + ... . It is clearly stated [2b, 3] 
that it is this constant X that must be zero, not the entire function ((rj), and not the con- 
stant C\ which has a different origin and requires reference to the growing solution. Martin 
and Schwarz [1] essentially repeat all this analysis and come to the same conclusion, but, 
regretably, portrait it as a Grishchuk's "mistake", rather than a complete agreement with 
his conclusion. After this clarification, we return to the "conserved" solutions determined 
by the constant C\. 

The authors of [1] do not analyze the physical meaning of ( and do not discuss its initial 
numerical value. They insist on the "nonemptiness" of the ( conservation law (see also |p~I[ ) 
and "show how to make use of the constancy of £" . They operate with the definition (|T2"D 
and arrive at their central statement (4.25): "During inflation w ~ — 1 and therefore the 
large amplification 

$„(t m ) ~ ^ \ ~, J _ 1 ji M U) (4.25) 
5 5 1 + w(ti) 

follows, which is the same as (4.7)". Formula (4.7) is 

$ m 22A + 3 1 2 1 

~ fa (4.7) 

$i 5 Pi + 1 1 + Wi 5 1 + Wi 

and is characterized as the "standard result" [1]. One more time Martin and Schwarz arrived 
at (4.7) in the form of their eq. (6.10). Apparently, the variant (6.10) is considered especially 
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valuable, because the authors of [1] note that they have now "obtained the 'standard result' 
(4.7) entirely within the synchronous gauge, without any reference to the constancy of £ 
or the joining conditions of Deruelle and Mukhanov". This last expression for the "large 
amplification" of $ directly goes over into the (incorrect) final conclusion: "the closer the 
inflationary epoch is to the de Sitter space-time, the less important are large-scale grav- 
itational waves in the CMBR today". So, we need to sort out this "large amplification" 
of $. 

Since we have in our hands a complete solution to the problem of evolution, it is not 
difficult to calculate the evolution of $ as well. The function $ is expressed in terms of the 
function Ji/a by flT3|). The dominant ("growing") term of Jt/a is a constant C\ (see (|Sp). 
So, in order to find the first nonvanishing approximation for we need to go to the next 
term in the long-wavelength decomposition of ft/ a. It was shown [3] that the more accurate 
solution for Ji/a reads 



a 



1 - n 2 ( ( aVr^dJ + C 2 [ ^- + ... (16) 
J GT7 J J a Ar y 



d?7 

a 2 7 



(this formula is given as (4.18) in [1]). Using (|i~6|) and neglecting the term with C2, we derive 



from (p~3|) the expression for $ which we will be working with: 



$ = -\cA 2 I a 2 jd V . (17) 



Let us first make a brief inspection of (|T7|) . It was shown [2] that the gauge-invariant 
potential $ is strictly zero at the de Sitter stage, that is, $ = for 7 = 0. This fact is reflected 
in the exact formula (|T3| ) and in the approximate solution fllT]) : the function $(77) = as 
long as 7(77) = 0. If 7(77) is not strictly zero but is very small, 7 -C 1, $(77) also is expected 
to be very small. If 7(77) is growing from zero, or almost zero, up to 7 = 2 or 7 = 3/2 (the 
cases of our interest) the function $(77) is also expected to grow. In other words, the absolute 
value of $(77) is expected to grow from zero or almost zero up to some finite number, which 
we will calculate shortly. The reheating transition can be either smooth or sharp, it does not 
matter. The integration of the possible step function 7(77) presents no difficulty, the function 
(|17|) remains continuous. This is simply a reflection of the exact result mentioned before: the 
function v, eq. (|TT|), and therefore (because a is continuous) the function $, eq. (P"3"|), must 
be continuous across a sharp transition. As expected, the first time derivative of $ must 
experience a finite jump at the sharp transition. This follows from the exact relationship 
$' = — a(7 + 1)$ — (l/2)a!7(/Z/a) which is a consequence of (|T3|) and (Q). Since the functions 
Ji/a and $ are continuous, the finite jump of <£>' is determined by the finite jump of 7. This 
finite jump of $' guarantees that even if the initial period of expansion was a strict de 
Sitter, the function $(77) would start slowly growing after the completion of that period, 
without developing any infinities, divergencies, or violations of linear perturbation theory, 
associated with the fact that 7 = at the initial stage. Correspondingly, the amplitudes 
in the Harrison- Zeldovich spectrum do not blow up to infinity. It was already emphasized 
[2] that there is nothing spectacular about the de Sitter case, the perturbation amplitudes 
remain finite and small. 

Let us now calculate (|PTD at the interval of expansion described by a power-law scale 
factor 
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<v) = Qv\ 1+P - 



(18) 



In this case, 7(77) is a constant: 7 = (2 + /3)/(l + (3) and 7 = ((3 — —2) corresponds to the 
de Sitter solution. The function (|i~7|) reduces to 



-C 

2 l7 3 + 2/T 



(19) 



Strictly speaking, the integration constant in (|17D produces an additional term proportional 
to a /a 2 , but this term is decaying and can be neglected, like the one (already neglected) 
which is associated with the constant C2 in flip]). The approximate solution (|19"D follows also 
directly from (|i~3"D and the exact solution for fi/a which can be found in the power-law cases 
©[2]. 

To get a complete and concrete result, we will now evolve the function $ through two 
power-law stages, and from one Hubble radius crossing to another. The first stage is 



the second stage is 



(20) 



The constants l m and rj e are so chosen that a (rj) and a{rf) are continuous at the transition 
point T) — t\\. The cases of particular interest are 7^ <C 1 and 7 TO = 2 or j m — 3/2. 
The function (Tl7f) takes the form 



3 + 2/3, 2 C ' lTiii * 



a Jr)i 



a 2 dr]. 



(21) 



The first term in (21): 



1 i + A 

"2^3+^ 



(22) 



gives the value of $(77) at the interval of evolution from the first crossing 77 = 7?j and up to 
t] = f]i. The second term in (^l|) gives the additional contribution at the interval of evolution 
from 77 = rji and up to the second crossing 77 = 7? m . If the first stage is close to the interval 
of the de Sitter expansion, fa m — 2, 7$ <C 1, then 



1 



4 



Ci(l + «;*). 



(23) 



Taking the integral (|2l|) with the help of 

$(?7m) 



1^ i+A 

-^7 



1 1 + 

3 + 2A 2 l7m 3 + 2/3,, 



we arrive at the final value of $(77): 

3+2/3„ 



1 



Jim ~ Ve 



Since a m (r) m ) > a m (?7i) and j m > 7,, the $(r? m ) simplifies: 



11 



2 x ""3 + 2/5 m 

Specifically for the matter-dominated era, f3 m — l,j m = 3/2, we get 



HVm) = ~Ci. (24) 



Thus, the gauge- invariant potential $(77) grows from its value (23) up to its value (0). The 
final value of $(77) is of the order of C\, i. e. $(77) barely reaches the numerical level of the 
characteristic metric amplitude ( [10| ) but never exceeds it. The perturbation 5p/po in the 
matter density will also reach only the level C\. 

The comparison of this calculation with the concept of the "large amplification" of scalar 
metric perturbations demonstrates that this concept is simply a misinterpretation. Imagine 
that we work with the function f(x) = 10 _20 sinx. Construct the ratio f(x m )/f(xi) where 
x m = vr/2 and Xi <C 1. The ratio is 1/xi, and it goes to infinity for X{ — > 0. It would be greatly 
misleading to think that the function f(x) had experienced a very "large amplification". The 
real meaning of this amplification is that the function f(x) can go through zero, but can 
never exceed the level 10 -20 . 

The situation with the "large amplification" of $(77) due to reheating [1] is similar. The 
exact formula ( |i~4"| ) for ( generates ( = —{1/2)C\ in the lowest nonvanishing approximation. 



The use of (|3j) and (H) shows that eq. (4.25) is trivially satisfied. The use of (gD, (H) 
and ( |23| ) shows that eq. (4.7) is also trivially satisfied. [To check that the intermidate part 
of (6.10) is trivially satisfied, one needs to use the continuity of the function ft/ a.] The 
particular way of evolution (as demonstrated above) of the gauge-invariant potential <E> is 
specific for this particular quantity. If, instead of $, we chose to work with a different 
potential, say 5> where $ = $/7, the idea of "large amplification" would not have even 
arisen. The quantity $ stays at the level C\ at the first crossing and does not surpass this 
level during all the time up to the second crossing. The potential $ is more adequate for the 
problem, because $ = (l/2n 2 )a(/2/a)', see eq. (|i~3|), and $ plays the role of the generalized 



"momentum" conjugate to the generalized "coordinate" Jija. 

Thus, the real meaning of the "large amplification" (4.7), (4.25), (6.10) is not the one 
that the authors [1] want to assign to it, with the associated incredible (incorrect) prediction 
of the larger and larger todays's amplitude of density perturbations for smaller and smaller 
1 + W{. The real meaning of this "large amplification" is the trivial fact that the function 
$(77) never goes above the level C\ (the constant level of the scalar metric perturbation 
h(r]), see eq. ([10])) even if $(77) has started from exceptionally small values at the almost 
de Sitter phase w,i ~ — 1. [It seems that the authors [1] agree with this fact at least in one 
part of the paper, in the discussion at the end of Sec. IVA, but their final conclusion is in 
conflict with this discussion.] One of inflationary misuses is to write the ratio of today's 
amplitudes for gravitational waves and density perturbations, in such a manner as if it was 
the contribution of gravitational waves which turned out to be, for some reason, exceptionally 
small in the limit W{ ~ —1 (see, for example, (7.1) in [1]). But the constancy of amplitude 
for gravitational waves has never been a matter of dispute. The reason for this small ratio 
is the claim about "big amplification" of the metric amplitude for density perturbations, on 
the route from the first crossing to the second crossing, which we have demonstrated to be 
incorrect. 
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After having sorted out all the questions of the evolution we need to return to the initial 
conditions which prescribe a numerical value to the constant C\. (The importance of initial 
conditions was emphasized in a recent work of Unruh ||16|| .) We need to say a few words 
about quantum normalization. 

The interest to primordial cosmological perturbations is related to the possibility of 
their generation from the inevitable zero-point quantum fluctuations by the mechanism of 
superadiabatic (parametric) amplification. This was first explored for the case of cosmolog- 
ical gravitational waves flOfl . The very first concrete example, a(t) oc t,a(n) oc e v , studied 
TUfl belongs to the class of background solutions which were later named inflationary. 



in 



There is of course nothing specifically inflationary in the mechanism itself. The generated 
amplitudes and spectral slopes depend on the evolution of the very early Universe and can 
be used for inferences about its behaviour. Specifically for the power-law scale factors ( |T8D . 
and in the long-wavelength approximation, the spectral characteristic amplitude h(n) (in 
logarithmic frequency intervals) was found to be h(n) « (lpi/l )n 2+ ^ . The generated field 
exists today in the so-called squeezed vacuum quantum state and carries specific statistical 
properties, but these properties are outside of the scope of our present discussion. The 
quantum considerations are only relevant to our discussion as much as they determine the 
constant C\ for gravity-wave solutions: C\ ~ (Ipi/Xi) (and Aj ~ l a for models with the 
parameter /3 not far away from the de Sitter value (3 = —2 ). 

Quantization of cosmological density perturbations is more complicated than quantiza- 
tion of gravitational waves, because the matter field fluctuations participate at the equal 
footing with the metric (gravitational field) fluctuations. Quantization of matter perturba- 
tions alone (for example, quantization of the scalar field perturbations in models considered 
here) would be as inconsistent quantization of the magnetic part of the electromag- 

netic field without quantization of its electric part. The model of a free quantized scalar 
field is usefull to the extent that it allows us to understand the way in which the fundametal 
constants G, c, h should participate in the final expression [2], but this model is certainly 
not sufficient for evaluation of the constant C\ for scalar metric perturbations. After all, it 
is the metric (curvature) perturbations that we are mostly interested in. We have to deal 
with a combined degree of freedom, which we call the "fundamental" field [2]. A spatial 
Fourier component of this field is essentially the dimensionless function fi/a, and its value 
(fi/a)i should be evaluated from quantum considerations. For models with 7 = const, that 
were actually considered [2] , the constant 1 / y/j combines with the constant in front of the 
relevant Bessel solution for fi(r]), as well as with the value cij of the scale factor a(n) at the 
first crossing for a given mode, to produce the dimensionless constant of our interest C\ for 
density perturbations. The order of magnitude normalization, based on the assigning of a 
half of the quantum to each mode of the "fundamental" field, gives C\ ~ (lpi/Xi) in full 
analogy with gravitational waves. And the spectral characteristic amplitude h(n) for the 
(scalar) metric perturbations is again h{n) ~ (lpi/l )n 2+l3 [2, 3]. [More accurate calculation 
with exact solutions gives some numerical preference to gravitational waves [2, 3].] 

If one does not want to deal with the "fundamental" field and prefers to quantize the 
(residual-gauge-invariant) scalar field perturbations which include metric perturbations, one 
is free to do that. The function that will do the job is ipi = tp\j\J\ + w + (3/4/t) 1 / 2 /i. 
Because of the exact relationship [2] </ , i( 7 ?) — (2/«)~ 1 ' 2 [/i/a — v/7/1], this function is ipi — 
(3/4/t) 1 / 2 (/2/a). On the other hand, the function ip 1 is the gauge- invariant scalar field variable 



13 



(called for example in the recent paper by (ps$) caculated in the synchronous "gauge". 



It is seen from these expressions that it would be incorrect to naively prescribe the quantum 
normalization to ipx on the grounds of the free scalar field quantization: (pi ~ Hi. If one 
takes this evaluation of the "vacuum fluctuations" and then sends 1 + w to zero, the value 
of , as well as the metric perturbation hi , become arbitrarily large already at the first 
crossing in proportion to + W{ . 

It is necessary to mention that sometimes one tries to make sense of the "standard" 
formula for (5p/p)f by saying that although there is not any "big amplification" on the way 
from one crossing to another, the formula itself applies only to cases ((p)i 3> Hi, that is, to 
\/l + Wi ^> 1 and yfy. ^> 1. If so, all the claims about dominance of density perturbations 
over gravitational waves are false from the very beginning, because they all are based on 
the opposite inequalities. In addition, this position seems to be a refusal to answer the 
question on what will happen in models with the opposite inequalities. The present author 
believes that there is no problem with quantization in those models as well, and that they 
were correctly treated in [2, 3]. In any case, the idea of a "big amplification" [1] in course 
of evolution, as an explanation for the "standard" formula, is incorrect and does not help in 
resolution of a possible new dispute - this time, over quantization. 

The conclusions of [2, 3] are relevant not only for the interpretation of the observed 
CMB anisotropies, but also, and may be more importantly, for the prospects of detecting 



relic gravitational waves by terrestrial and cosmic instruments [15 
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